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COMPARATIVE PLANETARY ATMOSPHERES: MODELS OF TRES-1 AND HD209458B 
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ABSTRACT 

We present new self-consistent atmosphere models for transiting planets TrES-1 and HD209458b. The 
planets were recently observed with the Spitzer Space Telescope in bands centered on 4.5 and 8.0 /im, for 
TrES-1, and 24 /mi, for HD209458b. We find that standard solar metallicity models fit the observations 
for HD209458b. For TrES-1, which has an T cS ~ 300 K cooler, we find that models with a metallicity 
3-5 times enhanced over solar abundances can match the la error bar at 4.5 /xm and 2a at 8.0/im. 
Models with solar abundances that included energy deposition into the stratosphere give fluxes that fall 
within the 2a error bars in both bands. The best-fit models for both planets assume that reradiation of 
absorbed stellar flux occurs over the entire planet. For all models of both planets we predict planet/star 
flux ratios in other Spitzer bandpasses. 

Subject headings: planetary systems, radiative transfer, binaries:eclipsing, stars: individual (TrES-1), 
stars: individual (HD209458) 



1. INTRODUCTION 

The detection of photons emitted by an extrasolar 
planet is a landmark feat in the history of astronomy. Ten 
years after the initial discovery of the first extrasolar giant 
planet (EGP) 51 Peg b (Mayor & Queloz 1995), Charbon- 
neau et al. (2005) report a dectection of infrared flux from 
TrES-1 and Deming et al. (2005b) report a similar detec- 
tion for HD209458b. Both observations were made with 
the Spitzer Space Telescope as each planet passed behind 
its parent star. While a number of papers have aimed 
at predicting the optical and infrared spectra of hot EGPs 
(Seager & Sasselov 1998; Marley et al. 1999; Barman et al. 
2001; Sudarsky et al. 2003), these results are the first in- 
frared detections (rather than upper limits) that can be 
compared with models. 

TrES-1, discovered by Alonso et al. (2004), is simi- 
lar to the well-known planet HD209458b (Charbonneau 
et al. 2000; Henry et al. 2000) in mass and orbital period 
but it receives less stellar irradiation and has a signifi- 
cantly smaller radius. Sozzetti et al. (2004) and Laughlin 
et al. (2005) have recently placed tighter constraints on 
the planet's mass and radius by improving our knowledge 
of the physical parameters of the parent star. 

TrES-1 and HD209458b are members of the class of "Pe- 
gasi planets" (or "hot Jupiters") that orbit their parent 
stars at less than ~ 0.1 AU and are presumed to be tidally 
locked. Unlike brown dwarfs of similar effective temper- 
ature, the luminosity of Pegasi planets is dominated by 
absorbed and re-radiated stellar flux, rather than intrinsic 
luminosity due to cooling of its interior. Atmospheric dy- 
namics will redistribute this absorbed energy around the 
planet with an efficiency that is currently unknown (Show- 
man & Guillot 2002; Cho et al. 2003; Burkert et al. 2005; 
Cooper & Showman 2005). Because of the dominance of 
absorbed flux over intrinsic flux, the planet's effective tem- 
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perature T c g , is equal to its equilibrium temperature, T eq , 
given by: 

1% = /(l - A)L+/(16irad 2 ), (1) 

where / is 1 if the absorbed radiation is able to be radiated 
away over the entire planet (4-7T steradians) or 2 if it only 
radiates on the day side (2ir steradians). A is the planet's 
Bond albedo, L± is the luminosity of the star, a is the 
Stefan-Boltzmann constant, and d is the planet's orbital 
distance. 

2. REVIEW OF THE OBSERVATIONS 

Charbonneau et al. (2005) used the Infrared Array Cam- 
era (IRAC) aboard Spitzer to observe the emission of the 
combined TrES-1 planet/star system. The timing of the 
observations included time before, during, and after the 
secondary eclipse (when the planet's light is blocked by 
the star). The observations were performed simultane- 
ously in bands 2 and 4, centered on 4.5 and 8.0 /xm, re- 
spectively. The planet+star system was indeed brighter in 
both bands when the planet was visible, indicating ther- 
mal emission from the planet was detected. The observed 
eclipse depths, in units of relative flux, were 0.00066 ± 
0.00013 at 4.5 fim and 0.00225 ± 0.00036 at 8.0 /xm. The 
HD209458b observations by Deming et al. (2005b) utilized 
this same method, except that their observation used the 
Multiband Imaging Photometer for Spitzer (MIPS) instru- 
ment at 24 /xm. The observed flux ratio was 0.00260 ± 
0.00046. 

From their observations Charbonneau et al. (2005) de- 
rived brightness temperatures in each of the two bands, 
T 4 . 5 = 1010 ± 60 K and T 8 . = 1230 ± 110 K. It is impor- 
tant to notice that Ts.o > T4.5, which was not predicted by 
models of Pegasi planets cited above, nor for brown dwarfs 
of similar T cS (Marley et al. 1996; Burrows et al. 1997; Al- 
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lard et al. 2001; Saumon et al. 2003a), nor is it observed 
in the mid infrared spectra of brown dwarfs to date (Roel- 
lig et al. 2004; Patten et al. 2004). Charbonneau et al. 
(2005) also derive a planetary T c ff = 1060 ± 50 K and 
A = 0.31 ± 0.14, but as the planet is not a blackbody, 
the values found are only suggestive. For HD209458b, 
Doming et al. (2005b) determine a brightness temperature 
T 24 = 1130 ± 150 K. 

3. METHODS 

To obtain our atmospheric pressure-temperature {P-T) 
profiles and spectra for the planets we employ a ID model 
atmosphere code that has been used for a variety of plan- 
etary and substellar objects. The code was first used to 
generate profiles and spectra for Titan by McKay et al. 
(1989). It was significantly revised to model the atmo- 
spheres of brown dwarfs (Marley et al. 1996; Burrows et al. 
1997; Marley et al. 2002), Uranus (Marley & McKay 1999), 
and EGPs (Marley 1998). It explicitly includes both inci- 
dent radiation from the parent star and thermal radiation 
from the planet's atmosphere. The basic radiative trans- 
fer solving scheme was developed by Toon et al. (1989). 
We use the elemental abundance data of Lodders (2003) 
and compute chemical equilibrium compositions following 
Fegley & Lodders (1994) and Lodders & Fegley (2002). 
In addition we maintain a large and constantly updated 
opacity database. We predict all cloud properties using 
the model of Ackerman & Marley (2001) with a sedimen- 
tation efficiency parameter f sec i = 3, which fits spectral 
observations of cloudy L-dwarfs (Marley et al. 2002). This 
model places 90% of the optical depth of a cloud within 
1 scale height of the cloud base. Further details can be 
found in Marley et al. (2002) and Marley et al. (in prep). 

We model the impinging stellar flux from 0.26 to 10.0 
/Ltm and the emitted thermal flux from 0.26 to 325 /im. All 
the relevant planetary parameters for TrES-1 are taken 
from Sozzetti et al. (2004). Those for HD209458b are 
taken from Brown et al. (2001). For the TrES-1 stel- 
lar model we use the Kurucz K0V model atmosphere 
computed for the Charbonneau et al. (2005) paper. For 
HD209458b we use a Kurucz (1993) GOV model (L* = 
1.6L Q ) with parameters described in Brown et al. (2001). 
The planet's radiative-convective P-T profile is arrived at 
iteratively until the net flux is conserved to at least one 
part in 10 6 . We compute all band-averaged flux density 
ratios using Eq. (1) from Charbonneau et al. (2005). Our 
model atmosphere code computes the P- T profile and low 
resolution spectra covering our full wavelength range. To 
generate a high-resolution spectrum we take the generated 
P-T profile and use a full line-by-line radiative transfer 
code, using the same chemistry and opacity database used 
in determining the P-T profile and low-resolution spec- 
trum (Saumon et al. 2000). 

4. RESULTS 

4.1. Atmospheric Pressure- Temperature Profiles 

Our computed atmosphere profiles arc shown in 
Figure 1. Two cloudless profiles for TrES-1 are shown, 
one under the assumption that absorbed stellar radiation 
is reradiated from the entire planet, and one for the as- 
sumption that the planet can only reradiate this energy 
on the day side. For these cases we derive a T e s of 1134 



K and 1352 K, respectively. Both of these profiles assume 
solar metallicity and an intrinsic temperature, T; nt of 250 
K. This Ti nt is the T e s the planet would have if it were 
in isolation. For the An case we also plot a profile for 
Tint = 100 K, which leads to a deeper radiative zone, but 
has no effect on the emitted spectrum. The presence of 
an outer radiative zone is predicted for the atmospheres 
of all Pegasi planets and is a caused by the dominance of 
external (stellar) radiation over internal heat flux (Guillot 
et al. 1996). The planet's current Ti nt can only be derived 
from thermal evolution models of the planet, which we do 
not compute here. For the An TrES-1 profile we show the 
model's range of brightness temperatures in 4.5 and 8.0 
fim bands. 

We also plot our derived P-T profile for HD209458b. 
This profile assumes An reradiation, solar metallicity, and 
T in t = 250 K. For HD209458b we find T eff = 1442 K, ~ 300 
K hotter than that for TrES-1. The model's brightness 
temperature at 24 /im is shown as a plus on this profile. We 
also plot the boundary where CO and CH4 have the same 
abundance. We predict CO will be the dominant (but not 
exclusive) carrier of carbon for both planets. Also shown 
are the condensation curves for Na2S, enstatite (MgSiOa), 
forsterite (Mg 2 Si04) and iron (Fe). Clouds form deep in 
the atmosphere of TrES-1 and therefore have no effect on 
the planet's spectrum. For HD209458b, we find that a po- 
tential enstatite cloud layer does affect the P- T profile and 
spectrum, as shown by the dash-dot profile. High clouds 
in the atmosphere of HD209458b are compatible with the 
available data from atmospheric transmission spectroscopy 
(Charbonneau et al. 2002; Fortney et al. 2003; Doming 
et al. 2005a), although the clouds in those studies are at 
mbar pressures. Our profile for HD209458b is very sim- 
ilar to the An model of Iro et al. (2004), although our 
profiles show hints of stratospheric temperature inversions 
due to absorption by near infrared water bands, which to 
our knowledge other authors have not found. 

4.2. Flux Ratios and Brightness Temperatures 

Our computed flux ratios for the HD209458b model at- 
mospheres match the Deming et al. (2005b) observation 
well. We find AF = 0.00305 for a cloudless model or 
with clouds of iron and forsterite, which is within the ler 
error bar. If the silicate cloud is not forsterite, but en- 
statite, which forms higher in the atmosphere, we find 
AF — 0.00316. Models computed assuming 2n reradia- 
tion would be brighter at 24 /im, leading to a poorer fit to 
the datum. 

For TrES-1 we will examine our model atmospheres in 
detail. In Figure 2 we show three panels of results. Panel 
A shows the surface flux of the model K0V stellar spec- 
trum along with two computed surface fluxes of planet 
TrES-1. These spectra are for 27r and An reradiation and 
solar metallicity. As expected, the 2n spectrum is brighter. 
Note that between 4-10 /im we find that emitted planetary 
flux is ~ 7-9 orders of magnitude brighter than reflected 
stellar flux. 

The lower two panels of Figure 2 show the planet /star 
flux ratio, as would be seen by a distant observer, which 
folds in the factor of ~ 60 difference in surface area of 
the planet and star. These panels compare our computed 
ratios to those measured by Charbonneau et al. (2005). 
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For 2tt reradiation our ratio at 4.5 /im is 4.5<r too high, 
indicating the planet model is far too luminous at these 
wavelengths. However, the fit at 8.0 /im is quite good. For 
the 47T model, the flux at 4.5 /im is a reasonable fit, falling 
just above la error bars, while at 8.0 /im the model falls 
below the 2a error bar. Both of these models have approxi- 
mately the same infrared spectral slope, which is obviously 
not as steep as Charbonneau et al. (2005) observed. This 
slope is mainly controlled by H 2 absorption on the blue 
side of the 4.5 /im band and strong CO absorption on the 
red side. At 8.0 /im, H 2 and CH4 are the dominant ab- 
sorbers. Based on these results, we rule out the 2tt model 
and examine the sensitivity of the spectral slope to model 
parameters. Our computed brightness temperatures for 
the 4?r models are T 4 . 5 = 1075 K and T 8 . = 995 K, while 
Charbonneau et al. (2005) found T 8 . > T4.5. Compared 
to a blackbody, the planet is redder while the models are 
bluer. 

A flux deficit at 4.5 /im, relative to atmosphere mod- 
els using equilibrium chemistry, is a common feature of 
Jupiter and T dwarfs (Golimowski et al. 2004). This is 
likely due to dredging of CO from deeper layers, which ab- 
sorbs flux from 4.5-5 /im (Fegley & Lodders 1996; Saumon 
et al. 2003b). We find this process cannot affect the atmo- 
sphere of TrES-1 as CO is already the dominant carbon- 
bearing molecule in the planet's atmosphere. 

We also varied the metallicity of the planet's atmo- 
sphere. We computed chemistry and opacity grids at 
metallicities of 3.2 QM/H]=0.5) and 5.0 ([M/H]=0.7) 
times solar. Recall that Jupiter's atmosphere is enhanced 
in heavy elements by factors of ~ 2-4. With 3.2 x or 5x 
solar metallicity the models fall within the la error bar 
at 4.5 /im and 2a at 8.0 /im. Panel C of Figure 2 shows 
the flux ratios for the 5 x solar abundance model. Increas- 
ing the metallicity by 3.2 x decreases the CH4/CO ratio by 
an order of magnitude. (The position of CH4/CO=l curve 
varies inversely proportional to the metal/hydrogen ratio.) 
This in turn increases the 8.0 /im band flux by weakening 
the strong CH4 band at 7.8 /im. Increasing the metallicity 
also increases the amount of H 2 0, which decreases the flux 
on the blue side of the 4.5 /im band. 

We also considered whether we could further improve 
the fit by increasing the C/O ratio from 0.5 (solar) to 0.7, 
therefore increasing CO absorption at 4.5 /im and decreas- 
ing H 2 absorption at 8.0 /im. However, increasing the 
C/O ratio (by increasing the carbon abundance) increases 
the CH4 abundance relative to CO (and H 2 0) and flux 
actually decreases in the 8.0 /im band. Increasing [M/H] 
from to +0.5 or increasing the C/O ratio from 0.5 to 0.7 
has very little effect on the strength of the CO band. This 
is because the band is formed very high in the atmosphere 
and is already saturated at solar abundances. 

Another avenue we pursued was increasing the tempera- 
ture of the upper atmosphere. In our models the pressures 
at which optical depth unity occurs in the 8.0 /im band are 
slightly smaller than those in the 4.5 /im band, although 
there is considerable overlap. (See Figure 1.) Therefore, 
increasing the temperatures at lower pressures preferen- 
tially increases the Tg.o relative to T4.5. We deposited 
additional energy in the form of a Chapman function (e.g. 
Chamberlain & Hunten 1987), into the upper atmosphere 
at pressures ranging from 1 to 10 mbar. Possible unmod- 



eled energy sources include conduction downward from a 
hot thermosphere (Yelle 2004), the breaking of gravity 
waves, or an unmodeled optical absorber, such as pho- 
tochemical products. Profiles that warm the atmosphere 
at 10 mbar by ~ 200 K (requiring ~ 10 7 erg cm~ 2 sec -1 , 
or 2.5% of the incoming stellar flux) generate more flux 
in the 8.0/im band, and reach the 2a error bars in both 
bands. Reprentative flux ratios of one of these models is 
shown in Panel C of Figure 2. 

In another trial (not shown) we modeled a strong con- 
tinuum opacity source from 2-6 /im. An appropriately 
strong absorber cuts down the flux in the 4.5 /im band 
while increasing flux in the 8.0 /im band, which allows the 
model to fit the la error bars in both bands. It is unclear 
what this mystery opacity source could be. Liang et al. 
(2004) found that HD209458b should be free of hydrocar- 
bon hazes as atmospheric temperatures are too high for 
hydrocarbons to condense. Their conclusion should hold 
for TrES-1, although it is ~ 300 K cooler. We note that 
the photochemical destruction of CH4 (Liang et al. 2003) 
would lead to slightly weaker absorption at 8.0 /im. 

5. DISCUSSION 

We find that a solar metallicity atmosphere for TrES- 
1 marginally fits the Spitzer observations. Models with 
enhanced metallicity are promising as they fall within the 
la error bar in the 4.5 /im band and 2a error bar in the 
8.0 /im band. Solar metallicity models with extra energy 
deposition into the stratosphere fit at the 2a level. As 
Table 1 shows, observations with IRAC at 3.6 /im look 
most promising for further separating the predictions of 
standard solar abundance models from others. For the 
model phase space we explore it is extremely difficult to 
obtain brightness temperatures of T 8 . > T4.5. One possi- 
ble remedy would be for the planet to have a cooler tro- 
posphere than we predict here, but a much hotter strato- 
sphere, which could cause the planet's absorption features 
to become emission features, and vice versa. We can re- 
cover this behavior with an ad-hoc P-T profile. We note 
that in Jupiter the 7.8 /im CH4 feature is seen in emission, 
not absorption. 

As shown above, standard solar abundance models fit 
the Doming et al. (2005b) MIPS observation well. How- 
ever, with this one data point metal enhanced atmo- 
spheres are not excluded. The known differences between 
HD209458b and TrES-1 are the former's larger radius and 
higher T c g. It is not yet clear if these atmospheres are 
substantially different in character. Possibilities include 
opacity due to silicate condensates in HD209458b but not 
TrES-1, and opacity due to photochemical products in 
TrES-1 but not HD209458b. Spitzer observations in ad- 
ditional bandpasses for both planets will help to charac- 
terize these atmospheres. Observations of the depth of 
the primary eclipse with Spitzer will allow the measure- 
ments of nightsidc fluxes of the planet, which would place 
strong constraints on atmospheric dynamics and on the 
energy redistribution between the day/night hemispheres. 
In Table 1 we collect our computed planet / star flux ratios 
for a variety of TrES-1 and HD209458b models in all four 
IRAC bands and the MIPS 24 /im band. 

Since models of the mid infrared spectra of brown dwarfs 
match observations to date (Roellig et al. 2004; Patten 
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et al. 2004), we suggest that while atmospheric composi- 
tion plays a role, a complete explanation for the red spec- 
tral slope of TrES-1 likely will involve energy from the 
parent star. Whether this is due to unknown photochemi- 
cal processes, unmodeled energy sources that dramatically 
modify the P-T profile, or some other process is not yet 
clear. 
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Fig. 1. — Pressure-temperature (P-T) profiles for TrES-1 and HD209458b. The solid lines are two profiles for TrES-1 assuming either 2-rr 
(hotter) or 4w (cooler) rcradiation of absorbed stellar flux. Both profiles assume T; nt = 250 K. For the An case, we also plot a profile with 
Tint = 100 K. For the 4n profile the thick black portion of the profile shows the extent of the brightness temperatures in the 8.0 fim band. 
The even thicker gray portion of the profile shows the brightness temperatures in the 4.5 (im band. We plot two 4w profiles for HD209458b. 
The dashed line is cloud-free and the dash-dot line includes the opacity of MgSiOa and Fe clouds. The model brightness temperatures at 24 
(im are marked with plusses. The thick dotted line is the boundary where CO and CH4 have the same abundance. Condensation curves for 
Na2S, MgSi03, Mg2Si04, and Fe are thin dotted lines. 



Table 1 



Model Flux Density Ratios in Spitzer Bands for TrES-1 and HD209458b 



Planet 


Model 


Toff 


3.6/im 


4.5^m 


5.8/xm 


8%m 


24//m 


TrES-1 


47T 


1134 


0.46 


0.81 


1.11 


1.48 


2.93 


TrES-1 


2tt 


1352 


1.06 


1.30 


1.69 


2.24 


3.69 


TrES-1 


47r,3.2x solar 


1144 


0.58 


0.78 


1.15 


1.62 


3.02 


TrES-1 


47r,5x solar 


1148 


0.61 


0.78 


1.17 


1.66 


3.06 


TrES-1 


47T, heating 


1165 


0.56 


0.90 


1.24 


1.71 


3.19 


HD209458b 


47T 


1442 


0.96 


1.12 


1.45 


1.90 


3.05 


HD209458b 


4-7T, cloudy 


1448 


1.02 


1.22 


1.56 


2.01 


3.16 



Note. — All ratios have been multiplied by 1000. Abundances are solar, unless 
noted. HD209458b "cloudy" includes enstatite and iron clouds. 
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Fig. 2. — Panel A: Emergent surface flux density (in erg s _1 cm -2 Hz -1 ) for the TrES-1 KOV parent star and two planet models — the 
higher curve assumes 2ir flux reradiation and the lower thick curve 4n. Relevant absorption features are labeled. Panel B: Planet/star flux 
density ratios for the models in panel A, as viewed by a distant observer. The full-width at half-maximum wavelength range of both bands 
are shown. The one- and two-cr error bars of the Charbonncau ct al. (2005) measurements arc shown as shaded boxes and empty boxes, 
respectively. The inset is a blowup of the 3-10 £tm region. Panel C: The gray curve is HD209458, the solid black curve is TrES-1 with 5 times 
solar metallicity, and the dash-dot curve is TrES-1 with solar metallicity and a heated upper atmosphere. (See text.) 



